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Abstract 



We introduce a family of correlated trial wave functions for the TV-particle 

ground state of an interacting Bose gas in a harmonic trap. For large N, the 

correlations lead to a relative energy decrease of a fraction compared to 

mean field Gross-Pitaevskii theory. The kinetic energy in the weakly confining 

direction turns out to be most sensitive to our correlations and, remarkably, 

is higher by as much as a few per cent for condensates with atom numbers of 

a few thousand. Thus, the predicted deviations from Gross-Pitaevskii theory 

originating from ground state correlations might be observed in momentum 

distribution measurements of small condensates. 
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Initiated by the first realizations of alkali Bose-Einstein condensates in magnetic traps 
a few years ago |1] there is now an ever growing number of experiments with atomic con- 
densates in laboratories all over the world. As a result of these efforts, a more and more 
detailed understanding of such quantum iV-body systems is emerging. Properties hitherto 
measured are successfully described by updated versions of many-body theories established 
many years ago, see || for a recent review. Very importantly, the 'condensate wave func- 
tion', the most relevant object well below the transition temperature, is determined by the 
mean field Gross-Pitaevskii (GP) equation |§. 

Among future challenges remains experimental access to the detailed nature of the N- 
particle state of these quantum systems well below the critical temperature, exploring new 
physics beyond the standard mean field description. This paper presents a theoretical backup 
of this aim, emphasizing non-mean-field effects in trapped condensates with a relatively 
small number of atoms. Consequences of a correlated ground state are expected to be 
most prominent in two-body or even higher-order correlation functions. Nevertheless, in 
this paper we show that correlations may already have a significant effect on experimentally 
directly accessible quantities like the kinetic energy of the condensate. This holds true, as 
we will show, whenever the particle number is not too large and the trap is anisotropic. 

The harmonic trap potential with frequencies u x = u y , and u z , provides the scales of 
the problem. We express energies in units of Ku) where u = (uj x u y u z ) s , the corresponding 
oscillator ground state length scale is d = y^j, and thus the momentum scale (h/d). In 
these units the Hamilton operator, representing the energy of N interacting atoms in an 
axially symmetric harmonic trap with anisotropy parameter A = u z /uj x , reads 

N ( n 2 1 "I 

Interactions between atoms are two-body collisions with a strength determined by the 
s— wave scattering length a of the atoms (measured in units of the oscillator length d). 
The problem thus depends on the three parameters N, A, and a. 

Bose-Einstein condensation occurs when the thermal de-Broglie wave length of the atoms 
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becomes larger than the mean distance between atoms. Then, loosely speaking, a notice- 
able fraction of all atoms occupy the same one-particle state. Near zero temperature, i.e. 
well below the critical temperature, we therefore expect the iV-particle state to be well 
approximated by a product state $ of one-particle states ip ||] , 

<&(r 1 ,...,r N ) =i/;(r 1 )i/;(r 2 )--^(r N ). (2) 

In second- quantized language with atom field operator ^(r), the product state (H) becomes 
a Fock state with particle number N, resulting in the more common expressions for the 
manifestation of a Bose condensate, for instance (ffl (r)ip(r')) = Ntp*(r)ip(r'). 

The total energy of the system, E = (ty\H\ty) with Hamiltonian ([!]), evaluated with the 
product state (|]) is the Gross-Pitaevskii energy functional 

E GP [i)} =N Jd 3 rr(r)[-^A + ^(X-I{x 2 + y 2 ) + \h 2 ) + 27ra(iV - l)|^(r)| 2 ]^(r). (3) 

The determination of the approximate iV-particle ground state is thus reduced to finding 
the minimum energy wave function of the energy functional resulting in the famous 
Gross-Pitaevskii equation || for the condensate wave function ip(r) (strictly speaking, in 
GP theory the factor (N — 1) in front of the interaction term is replaced by N; as we are 
interested in effects for a fixed, finite number of particles N, we stick to (N — 1) throughout 
this paper). 

By now it is experimentally established that the one-particle state ip(r) so obtained does 
in fact well describe the properties of the condensate. For these dilute bosonic gases at 
near-zero temperature, the assumption of a product wave function (0) for the iV-particle 
state is thus a good approximation. There is now also a rigorous proof available || showing 
that in the limit N — > oo, keeping the product Na fixed (dilute), the true iV-particle ground 
state energy is indeed given by the the minimum of Eqp[i/j]. 

Nevertheless, it is clear for finite N that due to two-body forces a product wave function 
(0) can only be an approximation for the iV-particle ground state and it is worth thinking 
about possible observable deviations for not too large atom numbers N. 



In particular, as the center-of-mass motion of the iV-body problem may be separated, we 
expect the true ground state to be a product wave function of center of mass r cm = J2i r i 
and relative coordinates, i.e. ^(r 1 , . . . ,r N ) = ?/> cm (r cm ) - ^(relative coordinates) with the 
wave function ip for the relative motion symmetric in the r$. This form of the wave function 
is clearly different from the product state (|2|). 

In this paper, we are not going to separate the center-of-mass motion, but instead we 
introduce new coordinates 

R^r. + ^-l)^, (z = l,...,iV) (4) 

that democratically incorporate the center-of-mass degree of freedom. The matrix C, with- 
out loss of generality, is chosen to be the diagonal matrix 

( n N 

C r (J 



C 



C r 



y c z j 



(5) 



where c r and c z are two positive parameters; C will play the role of an additional varia- 
tional variable and will help to lower the total energy. For a totally anisotropic trap one 
should introduce two parameters c x , c y instead of just c r . The Jacobian of transformation 
(|) is d^R = (detC)d Ar r, and notice that the choice C = 1 corresponds to the identity 
transformation Rj = r^. 

Apart from simplicity, transformation (ffl) is motivated by the fact that a wave func- 



tion ^(ri, . . . , r^v ) = VdetC $(Ri,...,Rjv) is a proper bosonic wave function whenever 
$(Ri, . . . , R^v) is a bosonic wave function in Rj coordinates. Crucially, the wave function 
obtained from a product wave function in Rj-coordinates 



tt(n, ...,r N ) = VdetC 0(R 1 )0(R 2 ) • • • <P(R N ) (6) 

will be a correlated bosonic wave function in atom coordinates unless C — 1, when (^) 
coincides with the usual product state (|]). 



4 



The actual values of the transformation parameters c r , c z and the shape of the wave 
function <p are fixed by the requirement that the total energy E = (^>\H\^>) evaluated in the 
space of correlated atomic wave functions @ should be minimal. It is straightforward to 
express the energy operator (|I|) of the N atoms in new coordinates Rj and the corresponding 
momenta. Variation with respect to C determines the parameters c r , c z of the transformation 
(H) with minimum energy, for a given wave function in (|6|). We find 

where Px = — an d, for instance, (Z 2 ) = ((f)\Z 2 \(f)). In deriving (|7|) we simplified using 
the fact that in determining the ground state \& we may restrict ourselves to wave functions 
0(R) with (R) = (P) = 0. 

The total energy E = fi?\H\ifj) based on the correlated state @ with the best matrix C 
from (0) becomes 

£ cor [0] = EarM - ~ (A-V(^ 2 + n- V^l + ^r)) 2 " \ {^^)~{^)) (8) 
which is manifestly lower than the uncorrelated Gross-Pitaevskii energy Eqp[<P] (^|) and one 
of the main results of this paper. 

Let us first discuss result (|8|) for the case of a very large number N of atoms. The 
difference between i? cor [0] and £"gp[0] is of the order of expectation values like (X 2 + Y 2 ) 
or (Pf), while the dominating term Pqp[0] is iV times larger. To leading order, therefore, 
we may evaluate the improved total energy with the state obtained from the minimum 
of Pqp[0] alone, i.e. with the solution of the usual GP equation. Since, for N — > oo, (P 2 ) 
('kinetic energy') may be neglected with respect to (R 2 ) ('potential energy'), we see from 
(|) that in this limit, the energy of the correlated state @ may be written as E COT [(j)} = 
^gp[0] — ^f-^gpM' where <fi is the solution of the usual GP equation and P^p the GP 
potential energy. Note that this energy decrease is possible due to the existence of the 
trap potential. Using the asymptotic Thomas-Fermi expressions for the total GP energy 
Pgp = 7j (15(JV - l)aY and the GP potential energy E%$[<j>] = (15(iV - l)o)» §, we 
find for the energy of our correlated approximate ground state (||) 
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£ C or = (1 - ^) E gp as N - oo, (9) 

a small decrease for the assumed large N. 

Far more relevant are effects of the correlated ground state (||) on the properties of small 
condensates. Let us therefore concentrate on the full expression ([§]) for the total energy 
E C ot[<P] of the correlated state. As in usual GP theory we may proceed analytically 

and use Gaussian trial wave functions 0(R) parameterized by two parameters S r and E 2 
for the radial and axial widths (note that these are widths in R coordinates and should not 
be confused with the condensate widths in physical space, the difference being of the order 
N^ 1 ). The full energy functional (§) becomes a function of S r and only, and it is a simple 
numerical task to find its minimizing values. Results based on this Gaussian approximation 
are shown in the following Figures. 

In Fig.l we show the parameters c r (Fig. la) and c z (Fig. lb) of our transformation on 
correlated coordinates (f|) for the minimum energy correlated state (|6]) as a function of parti- 
cle number N, and for two different anisotropies A = 0.04 (cigar), and A = 3 (pancake). The 
three corresponding different graphs in each figure represent different interaction strengths: 
a = 0.004 (solid line), a = 0.008 (dashed line), and a = 0.012 (dotted line). Recall that 
c r = c z = 1 would correspond to the usual uncorrelated GP ground state (H). We clearly see 
that the transformation parameters c r and c z are larger in the weakly confining direction, 
and become larger for larger interaction strength a. 

Is there any hope to measure effects induced by the correlations of the improved ground 
state @? Among the various contributions to the total energy of the condensate, it turns 
out that the small kinetic energy in the weakly confining direction is most sensitive to our 
correlations and exhibits an increase of a few per cent compared to GP theory, for small 
condensates. This increase is compensated for by a larger decrease of potential energy of 
the correlated ground state, such that the total energy is indeed lowered. 

In Fig. 2 we show the kinetic energy in the weakly confining z-direction E^ ,z of the 
correlated state (|j) as a function of atom number N for an anisotropy A = 0.04 (cigar). As 
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in Fig.l, results are shown for three different interaction strengths, a = 0.004 (solid line), 
a = 0.008 (dashed line), and a = 0.012 (dotted line). Recall that energies are measured in 
units of fuD, i.e. the kinetic energy shown in Fig. 2 is very small compared to the potential 
and internal energy of the condensate. 

In Fig. 3 we compare the predictions of correlated (ffl) and uncorrelated (Q) ground state 
for the kinetic energy E km ' z in the weakly confining ^-direction, again for a cigar shaped 
condensate with A = 0.04 (both quantities evaluated in a Gaussian approximation). Shown 
is the relative energy increase 

A £kin,z = ( £ km,z _ ^km.z^km^ compared tQ Qp theory . 

We see that the difference may be a few per cent, depending on interaction strength and 
atom number. For larger anisotropy, the observed difference is even larger. Similar, yet 
less pronounced results are obtained for the change of radial kinetic energy in pancake-type 
condensates. 

The difference being a few per cent, high precision momentum distribution measurements 
might be able to distinguish between the predictions of correlated and uncorrelated ground 
state. In principle, the kinetic energy could be measured indirectly via the virial relation 
£kin, z = E rci _ £ P ot,r from the knowledge of release energy E rel = E km + E int § 

and radial potential energy E pot,T . In this approach, however, the small kinetic energy is 
measured as the difference of two very large quantities and it seems doubtful whether the 
required precision can be achieved. 

Far more promising are direct methods to measure the momentum distribution of the 
condensate, as recently established through Bragg spectroscopy 0. The high precision 
achieved in this experiment might be sufficient to confirm the predicted deviations from GP 
theory experimentally. 

Let us summarize this paper. We found a correlated iV-particle wave function (|6]) for 
an interacting Bose gas in a harmonic trap with lower total energy than the mean field GP 
product state. The existence of the trap potential is crucial for the correlations; the energy 
decrease is roughly the potential energy of one atom. We have concentrated on effects on 
simple averaged quantities, like kinetic, potential or release energy. Consequences of our 
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correlations will be most significant for relatively small condensates and a more detailed 
investigation of further effects is required, focusing on higher-order correlation functions. 

Studying the dynamics of the condensate in a correlated state \& will also be of interest 
as excitation frequencies can be measured with fairly high precision. Based on a hydrody- 
namic approach to superffuids, non-mean field corrections to the frequency of elementary 
excitations due to a finite gas parameter have been calculated in [|lOj in the large iV-limit. 
Note that for dynamics, we have to replace the correlated energy functional by the more 
general expression, valid for nonvanishing (R), (P). 

The correlations reported in this paper depend on the finite number of atoms and the 
existence of the trap; both conditions are met by current experiments. They are thus of 
very different nature than zero-temperature non-mean field effects described by Bogoliubov 
theory (quantum depletion) @JTl|]. In our approach, Bogoliubov theory in i?-coordinates 
may be introduced on top of the correlated ground state (^). We conclude that there is no 
obvious relation between the correlations described by our wave function (|6]) and the usual 
Bogoliubov quantum depletion corrections. 

Last but not least, recent progress in high precision spectroscopy of the momentum 
distribution of the condensate || might lead to an experimental confirmation of our results. 

I would like to thank R. Graham, M. Fliesser, and J. Reidl for numerous fruitful discus- 
sions. Support by the Deutsche Forschungsgemeinschaft through SFB 237 "Unordnung und 
groBe Fluktuationen" is gratefully acknowledged. 
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Fig. 1. Parameters c r (Fig. la) and c z (Fig. lb) of our transformation on correlated co- 
ordinates for the minimum energy state $ as a function of particle number N for two different 
anisotropies A = 0.04 (cigar), and A = 3 (pancake). The three corresponding different graphs in 
each figure represent different interaction strengths: a = 0.004 (solid line), a = 0.008 (dashed line), 
and a = 0.012 (dotted line). The uncorrelated case is c r = c z = 1. 
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Fig. 2. Kinetic energy E^ ,z of the correlated state in the weakly confining z-direction as 
a function of atom number N for an anisotropy A = 0.04 (cigar). As in Fig.l we show results for 
three different interaction strengths, a = 0.004 (solid line), a = 0.008 (dashed line), and a = 0.012 
(dotted line). 



12 
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Fig. 3. Relative difference between the predictions for the kinetic energy E km,z in the weakly- 
confining z-direction based on the correlated state and the uncorrelated GP state as a function 
of atom number N for an anisotropy A = 0.04 (cigar). As in Fig.l and 2, we show results for 
three different interaction strengths, a = 0.004 (solid line), a = 0.008 (dashed line), and a = 0.012 
(dotted line). 
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